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New global periodic orbit coUision/separatrix reconnection scenarios in the standard nontwist 
map in different regions of parameter space are described in detail, including exact methods for 
determining reconnection thresholds that are implemented numerically. The results are compared 
to a break-up diagram of shearless invariant curves. The existence of meanders (invariant tori that 
are not graphs) is demonstrated numerically for both odd and even period reconnection for certain 
regions in parameter space, and some of the implications on transport are discussed. 



In recent years, area-preserving maps that 
violate the twist condition locally in phase space 
have been the object of interest in several studies 
in physics and mathematics. These nontwist maps 
show up in a variety of physical models, e.g., in 
magnetic field line models for reversed magnetic 
shear tokamaks. An important problem is the 
determination and understanding of the transi- 
tion to global chaos (global transport) in these 
models. Nontwist maps exhibit several different 
mechanisms: the break-up of invariant tori and 
separatrix reconnections. The latter may or may 
not lead to global transport depending on the 
region of parameter space. In this paper we 
conduct a detailed study of newly discovered 
reconnection scenarios in the standard nontwist 
map, investigating their location in parameter 
space and their impact on global transport. 



I. INTRODUCTION 

In this paper we consider the standard nontwist map 
(SNM) M, as introduced in Ref. 



Xn+l =Xn+a{l- yl^i) 
Vn+i = Vn - 6sin(27rx„) , 



(1) 



where (x,y) G T x R are phase space coordinates and 
a, 6 £ M parameters. This map is area-preserving and 
violates the twist condition. 



^0 y{xi,yi), 



(2) 



along a curve in phase space, called the nonmonotone 
curve. |3| Some basic concepts used throughout the paper 
are reviewed in Appendix IXI 

Nontwist maps are used to describe many physical sys- 
tems, e.g., magnetic field lines in tokamaks (see, e.g., 
Refs. nil mill 1113) and stellarators[ll[l3 Masma 
physics); planetary orbits. [T2| stellar pulsations |13j (as- 
tronomy); traveling waves. |llll4l| coherent structures and 



self-consistent transport |l5l| (fluid dynamics). Additional 
references can be found in Refs. [lil ll^ . Apart from 
their physical importance, nontwist maps are of mathe- 
matical interest because important theorems concerning 
area-preserving maps assume the twist condition, e.g., 
the KAM theorem and Aubry-Mather theory. Up to now 
only few mathematical results are known, [l&i Il9l l20l l2ll 
|23 | Recently, it has been shown[2l 0| that the twist 
condition is violated generically in area-preserving maps 
that have a tripling bifurcation of an elliptic fixed point. 
For studies of nontwist Hamiltonian flows, we refer the 
interested reader to Ref. and references therein. 

Remark 1. Although the SNM is not generic due to its 
symmetries (see Appendix , it captures the essential 
features of nontwist systems with a local, approximately 
quadratic extremum of the winding number profile. 

Nontwist maps of the annulus exhibit interesting bifur- 
cation phenomena: periodic orbit collision and separatrix 
reconnection. The former, which applies specifically to 
collisions of periodic orbits of the same period, such as 
the so-called up and down periodic orbits that occur in 
the SNM, can be used to calculate torus destruction. [T^ 
The latter is a global bifurcation when the invariant man- 
ifolds of two or more distinct hyperbolic orbits with the 
same rotation number connect leading to a change in the 
phase space topology in the nontwist region. We briefly 
review previous studies of reconnection in nontwist sys- 
tems. 

Howard and Hohs|23 deflned a quadratic nontwist 
map closely related to Eq. ^ and studied numerically 
the reconnection of low-order resonances exhibiting ho- 
moclinic and vortex-like structures. Defining an aver- 
age Hamiltonian they predicted the reconnection thresh- 
old for period-one and period- two fixed points. Howard 
and Humphcrys extended the study to cubic and quar- 
tic nontwist maps.|23 These reconnection scenarios had 
been conjectured by Stixd in the context of the evolu- 
tion of magnetic surfaces in the nonlinear double-tearing 
instability, and were seen by Gerasimov et a/.'28| in a 
two-dimensional model of the beam-beam interaction in 
a storage ring. 

The first systematic stud y o f reconnection was done 
by van der Weele et aL|24L l29l | in the context of area- 



2 



preserving maps with a quadratic cxtrcmum. As far as we 
know the terminology "nontwist" and "meanders" origi- 
nates from there. 

Del-Castillo-Negrete et aZ.^3l devised an approximate 
criterion for the reconnection threshold of higher-order 
resonances based on matching the slopes of the unstable 
manifolds of the reconnecting hyperbolic orbits. 

Continuing the work of Carvalho and Almeida. |30| 
Voyatzis et aL|3ll| studied reconnection phenomena in 
nontwist Hamiltonian systems under integrable perturba- 
tion with cubic winding number profiles. Applying Mel- 
nikov's method, they showed the transverse intersection 
of manifolds for arbitrarily small nonintegrable pertur- 
bations for one reconnection scenario when the winding 
number has a local extremum. |32j 

The influence of manifold reconnection on diffusion was 
studied in the r egio n of strong chaos by Corso et al. in a 
series of papers. [I^, 0, iMl 

A very different criterion for the reconnection thresh- 
old was proposed by Petrisor in Refs. ^E^: If two hy- 
perbolic orbits have a heteroclinic connection, their ac- 
tions coincide. As noted in Ref. |37j . for odd-period hy- 
perbolic orbits in the SNM this criterion reduces to the 
action being zero at the point of reconnection. This crite- 
rion was implemented numerically in Ref. to estimate 
some reconnection thresholds for odd-period orbits in the 
SNM. 

Remark 2. It has been noted that the above result about 
the equality of actions of reconnecting hyperbolic orbits is 
only approximately true in the ncar-integrable limit, [s^ 

For completion, we mention a few related studies: re- 
connection phenomena and transition to chaos in the 
Harper map,^3j ^| degenerate resonances in Hamil- 
tonian systems with 3/2 degrees of freedom, ^3 a-nd 
zero dispersion resonance in the study of undcrdamped 
oscillators, [isj 

Key to the analytical and numerical exploration of 
the standard nontwist map is the map's invariance un- 
der symmetries, reviewed in Appendix |b1 Of particular 
significance are the indicator voints.]4^ fixed points of 
some of the symmetries of the SNM, whose importance 
was first recognized by Shinohara and Aizawa. These 
points were independently re-discovered by Petrisor in 
the analysis of the reversing symmetry group |45| of non- 
twist standard- like area-preserving maps. In Ref. |46l |. 
it was shown that a shearless invariant torus crosses the 
X-axis at two points. This led the authors to devise a 
criterion to determine the approximate location in the 
(a, 6)-parametcr space of the break-up of shearless in- 
variant tori for many winding numbers (sec Sec. IIIEI) . 

Subsequently, based on numerical observations. Shino- 
hara and Aizawa used indicator points to propose ex- 
act expressions for the collision threshold of even-period 
orbits and a method to determine numerically the recon- 
nection threshold for odd-period hyperbolic orbits. 

The goal of the present paper is to describe 
reconnection-coUision phenomena in detail in all regions 
of (a, 6)-paramcter space for the standard nontwist map. 



The details of the two main scenarios (odd-period and 
even-period) depend crucially on the (a, 6)-region, and 
have not been discussed exhaustively up to now. The 
main tool we use here is the numerical implementation 
of two criteria for collision thresholds: the analytic one of 
Ref. 01 and the numerical one of Ref. ■ The result- 
ing curves are compared with the (a, 6)-space break-up 
diagram (see Sec. IIIEj) . produced by a modified version 
of Refs. [3j EBI ■ Early results of this investigation were 
reported in Ref. , and some of them will be elaborated 
below for completeness. 

The paper is organized as follows. We review some ba- 
sic concepts of nontwist systems and the SNM in SecUll 
Some novel reconnection and collision scenarios for orbits 
of even and odd periods are described in Sec. IIII Al and 
Sec. IIII Bl resp ectively, applying and extending methods 
from Refs. The results arc discussed in the con- 

text of the break-up diagram in Sec. II VI In Sec. Ivl we 
give our conclusions and indicate some directions of fu- 
ture research. The appendices contain basic definitions 
and a brief summary of symmetry properties of the SNM. 

II. REVIEW OF NONTWIST MAPS 

In this section, we review some fundamental concepts 
that are required for further explorations of the SNM. 
Fig. n] shows a typical phase space plot along with the 
plot of winding number lo versus y (henceforth called 
winding number profile) for the central section of the y- 
axis. 

The phase space of the SNM consists of two twist 
regions, "above" and "below" the center, i.e., regions 
in which Eq. Q is satisfied. Since the twist condi- 
tion is violated along the nonmonotone curve C given by 
y = 6sin(27ra;), only orbits with points falling on C, re- 
ferred to as nonmonotone orbits by Petrisor, and orbits 
with points on both sides of C, called pseudomonotone, 
arc affected by the nontwist property. Note that, for any 
6^0, the curve C is not an invariant torus. 

Like the phase space plot, the winding number profile 
in the outer regions looks the same as that of a twist 
map, showing, e.g., the familiar plateaus associated with 
islands around periodic orbits (here orbits with wind- 
ing numbers 3/5, 7/12, and 4/7). The only distinctly 
nontwist effect, aside from the existence of the overall 
maximum (here at 3/5) and of multiple periodic orbits 
of winding numbers less than the maximum, is the valley 
below the 3/5 plateau, which gives rise to a variety of 
phenomena discussed in Sec. IIIII 



A. Reversing symmetry group and shearless curve 

When the winding number profile has a local extremum 
at an irrational value of tj, the corresponding invariant 
torus is called shearless. As discussed by Petrisor, 2] the 
SNM has at most one homotopically nontrivial shearless 
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FIG. 1: Example of the phase space for standard nontwist 
map (top) and winding number profile along y-axis (bottom) 
at a = 0.615, b = 0.4. The symmetry lines, nonmonotone 
curve C, ^-invariant curve 7s, and indicator points are dis- 
played. 



invariant torus 75 that is also invariant under the revers- 
ing symmetry group G (reviewed in Appendix IB|) . When 
it exists, 7s also contains the indicator points 

4°'^^ = (44),and.r=(^Tio), (3) 

and all their iterates. As we will discuss in Sec. 11111 
in some regions of parameter space, there exist several 
shearless tori, but only one of them is ^-invariant. 

The significance of a shearless torus is that it acts as a 
barrier to transport, whereas the nonmonotone curve C 
does not. Fig. ^ shows the nonmonotone curve and the 
shearless torus, along with symmetry lines and indicator 
points. 

The standard nontwist map is non-generic because it is 
time-reversal symmetric as well as invariant under a sym- 
metry, as reviewed in Appendix^ Nevertheless, most of 
the phase space phenomena described in this paper are 
also observed in arbitrary nontwist systems, even though 
the exact definitions of many of the concepts introduced 
to study them cannot be generalized to these systems. 



B. Periodic orbit collision and standard separatrix 
reconnection 

As mentioned earlier, one consequence of the viola- 
tion of the twist condition is that the SNM has more 
than one orbit (either invariant tori or chains of peri- 
odic orbits) of the same winding number. These orbits 
can collide and annihilate at certain parameter values. 
The collision of periodic orbits involves another purely 
nontwist phenomenon, namely the reconnection of the 
invariant manifolds of the corresponding hyperbolic or- 
bits. These reconnection-collision sequences in the SNM 
are distinctly different for orbits of even and odd periods. 
Here we will only give a brief account of the simplest ver- 
sion of both sequences, which we will refer to as the stan- 
dard scenarios. |29| We show phase space plots and wind- 
ing number profiles for several steps in these sequences. 
The description of more intricate reconnection-collision 
scenarios appears m Secs. lIII Al and llll Bl Changing the 
parameter values of the SNM we see the following se- 
quences: 



• For orbits with even period, i.e., with the same sta- 
bility type for both the up and down orbits on a 
symmetry line (Fig.|2ti): collision of the hyperbolic 
orbits which is also the threshold for reconnection 
(Fig. 1213); the "dipole topology" in which the hy- 
perbolic orbits have moved off the symmetry lines 
(Fig. 12:); the collision of elliptic orbits that coin- 
cides with the collision of the non-symmetric hy- 
perbolic orbits (Fig. |3i) leading to annihilation of 
these periodic orbits (Fig. |2t). 
The winding number profile shows a maximum that 
is greater than, equal to, and less than the winding 
number of the periodic orbit before, during, and 
after this process, respectively. 



• For orbits with odd period, i.e., with opposite sta- 
bility type for the up and down orbits on a sym- 
metry line (Fig. ^Cj: reconnection of hyperbolic 
manifolds of up and down orbits (Fig. Ot); ap- 
pearance of non-KAM meandering orbits (elabo- 
rated on in Sec. lIIDll . homoclinic separatrices, and 
dimerized chains (Fig.|3;); hyperbolic-elhptic colli- 
sion (Fig. |3Jl) leading to annihilation of these peri- 
odic orbits (Fig.EJ;)- 

As above, the winding number profile shows a 
global maximum that is greater than, equal to, and 
less than the winding number of the periodic or- 
bit before, during, and after this process, respec- 
tively. But in addition, there is a local minimum 
- associated with the appearance of the meander- 
ing orbits - which persists even after the collision. 
(See SecHTDl) 
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C. Bifurcation and indicator curves 

The parameter values for the threshold of collision of 
periodic orbits with a fixed winding number are numeri- 
cally observed to lie on a smooth curve in (a, 6) -parameter 
space, called bifurcation curve, which was first defined in 
Ref. 

The m/n-bifurcation curve b = $m/n,i(o) is the set of 
(a, b) values for which the m/n up and down periodic or- 
bits on the symmetry line Si (see Appendix ^ are at the 
point of collision. (For clarity, the figures in this paper 
denote these curves by "bcl," "bc2" etc. for symmetry 
lines si, S2 etc.) The main property of this curve is that 
for (a, 6) values below b = $m/n,i(fl)i the r/s periodic 
orbits, with r/s < m/n, exist. Thus, m/n is the max- 
imum winding number for parameter values along the 
m/n-bifurcation curve. 

For n odd, bifurcation curves along si through S4 co- 
incide in the SNM because of the map's high degree of 
symmetry. For n even, bifurcation curves along si and S2 
are separate from the ones along S3 and S4 for any b 0. 
Shinohara and Aizawafi^ used the numerical observa- 
tion that at the point of hyperbolic collision (Fig.|2l3) for 
even-period orbits (but not for the odd-period case) . two 
of the indicator points belong to the hyperbolic periodic 



orbit. This implies that (for period 2n) 

^^2„^(0.1) ^ ^(0,1) 

for either j = or j = 1, where 2^°'^'' are given in Eq. 
The symmetries of the map further imply that the indi- 
cator points map onto each other after n iterations, i.e., 

M-zf^'^^z\'^'K (5) 

By solving these equations for the two unknowns {a,b), 
we can obtain exact expressions for the bifurcation 
thresholds. Some of the resulting curves for low-period 
orbits are given in Table U They are related to the ones 
given in Ref. |44j by a simple transformation. 

We will later see that when more than two chains of 
periodic orbits exist, the collision threshold for some of 
them, but not for the others, is given by the above crite- 
rion. Hence we will find it useful to introduce the notion 
of indicator curves b = '^m/nio,) defined by Q (or ||SJ)) . 
When only two chains of m/n-orbits exist, either ^m/n,i 
or $m/n,3 coincides with ^'m/„. Hence we label indicator 
curves by symmetry lines, e.g., ^'m/n,i- (We will denote 
these curves by "icl" and "ic3" in the figures for clarity.) 
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TABLE I: Some exact expressions for indicator curves for even-period orbits. 



Remark 3. Shinohara and Aizawa also proposed a numer- 
ical criterion to find the odd-period reconnection thresh- 
old. They discovered numerically that at the threshold 
successive iterates of indicator points approach the same 
hyperbolic periodic point of the reconnecting chains. For 
details sec Ref . . 

D. Meandering orbits 

Another characteristic of nontwist maps is the occur- 
rence of meandering orbits, which are readily observed 
in the standard reconnection-coUision scenario for odd- 
period orbits, but not in the one for even-period orbits. 
(For non-standard even-period reconnection-coUision sce- 
narios in which meanders occur, see Sec. IllliO ) As seen 
in Fig. 121:; (and also in Fig.^, in the region surrounding 
7s, confined by two dimcrized chains, new periodic orbits 
and non-KAM tori appear, i.e., orbits and tori that did 
not exist at zero perturbation (b = 0). These tori are not 
graphs over the x-axis and have been called meanders or 
meandering curves. |22l Such invariant tori can occur 
only in nontwist maps because any invariant torus for a 
twist map must be a graph over x. 

It is observed numerically that in the meandering re- 
gion the winding numbers of the meandering orbits are 
less than the winding number of the reconnecting peri- 
odic orbits. Conversely, however, a "valley" in the wind- 
ing number profile does not imply the existence of mean- 
ders as seen, e.g., in Fig. The existence of the local 
minimum and the valley leads to four or more chains 
of periodic orbits for certain winding numbers. We will 
discuss in Sec. IIII Al and Sec. IIII Bl thc reconnection and 
bifurcation phenomena in this scenario, which we call the 
non-standard scenario. 

E. Break-up diagram 

Since the shear less curve 75, whenever it exists, poses 
a barrier to transport between the two twist regions, 
studying its break-up is of considerable practical interest. 
Highly detailed studies of break-up of shearless curves 
were conducted for a few winding numbers using Greene's 
residue criterion. 0, 0, ^3 Though the method is very 
precise, it is not suitable for an exploration of all param- 
eter space. Shinohara and Aizawa[46j obtained a rough 
estimate for the break-up threshold of many shearless 



curves by investigating for a range of parameter values 
whether iterates of one of the indicator points remain 
bounded. 

Here, we implement the following slightly different 
strategy, using the fact that if the winding number of 
the orbit of any point exists, then the orbit is not chaotic 
- it is cither periodic or quasiperiodic. We calculate the 
sequence uii = Xi/i for the iterates {xi,yi) = A/P{z) of 

one of the indicator points z^j*''^'^ ■ The winding num- 
ber is assumed to exist if we can find some N such that 
\uji - uji+i\ < e, sup„<jv > i^i, and inf„<jv i^n < i^i 
foi N < i < N + M. We use z = (a/2 1/4,0), e = 
lO""^, M = 10■^ and the maximum N used is 2.9 x 10^. 
If the winding number sequence displays larger fluctua- 
tions, we assume the orbit to be chaotic, i.e., the torus 
to be destroyed. This criterion is also only approximate, 
since the value of the winding number might converge 
for the number of iterations used, but further iterations 
would reveal fluctuations, or vice versa. 

The method of Ref. and ours deliver similar results. 
However, it seems that the computation of the winding 
number provides better means of monitoring and con- 
trolling its accuracy (aside from giving us the winding 
number of the shearless curve as a useful side-product). 

The boundary of the resulting break-up diagram dis- 
plays a fractal-like structure (see, e.g.. Fig. Ill|l . The 
analysis of Refs. [Til ItR li^ indicates that the highest 
peaks correspond to the break-up of shearless invariant 
tori with noble winding numbers (see Refs. 0, H^l for 
discussion). We will comment in Sec. llVl on the relation 
between the break-up diagram, reconnection, and the in- 
dicator and bifurcation curves. 



III. NON-STANDARD RECONNECTION AND 
COLLISION SCENARIOS 

As noted in Sec. IllDl the winding number profile can 
have a valley and a local minimum (denoted by tiirnin) 
between two local maxima (denoted by Wmax < ^max)- 
(In many cases it is observed that t^max = "-^maxj ^^^^ 
not distinguish between the two maxima.) Thus there 
are four orbits for each winding number in the range 
^min < LO < Wmax- The maxima and the minimum are 
seen to decrease when the perturbation, b, is increased. 
This gives rise to the following scenarios for reconnection 
and collision of m/n-pcriodic orbits: 
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• For Wniin > m/n^ there are two chains of m/n- 
periodic orbits. 

• When Wmin reaches the rational value m/n^ peri- 
odic orbits of winding number m/n are born and 
subsequently reconnect, in addition to already ex- 
isting orbits. We call the new orbits inner orbits 
while the already existing ones will be called outer 
orbits. 



0.245 - ic3 

be 1 - 
be 3 



• When (jJmax reaches m/n, the inner orbits reconnect 
and collide with the outer orbits. 

• For Wmax < rnjn, no ?7i/7i-periodic orbits exist. 

The reconnections and collisions that occur in the above 
scenario are locally the same as those seen in Sec. IIIBI 
But because of the presence of four chains of periodic 
orbits, the global topology is considerably more compli- 
cated as will be seen below. 

Remark 4. A valley in the winding number profile shows 
up after reconnection, together with meandering orbits, 
and also persists after the collision of the reconnecting 
orbits (Fig. The scenario described in the above 

paragraph typically occurs for parameter values slightly 
above (in b vs. a parameter space plots) the bifurcation 
curves of odd-period orbits. 

Remark 5. When a valley in the winding number profile 
exists, the winding number of the shearless curve is the 
local minimum ujuiim if ^min is irrational. 

Remark 6. When tJmin = m/n, the reconnection process 
of the inner (meandering) orbits involves orbits, called 
second- order meanders.\2^ that meander around these 
inner orbits. The winding number profile shows a "hill" 
at the bottom of the valley. Such a process might give 
rise to arbitrarily higher order meanders. |2, 123 . l38| 



A. Non-standard scenarios for even-period orbits 

A first hint that for the non-generic SNM the even- 
period reconnection scenario can be more complicated is 
found in Ref . [4^ . In some region of (a, 6)-parameter 
space, phase space portraits show the appearance of me- 
anders and additional periodic orbits. For most regions of 
(a, 6)-space the indicator and bifurcation curves are seen 
to coincide. But in regions of parameter space where the 
winding number profile has a valley, bifurcation curves 
and indicator curves separate (and can even cross each 
other) leading to various reconnection scenarios and the 
appearance of meanders. 

This is seen most clearly near a = 1, which is the bi- 
furcation curve for the 1/1-periodic orbits, but, as men- 
tioned above, occurs above every odd-period bifurcation 
curve. We discuss in detail the example of 7/8-periodic 
orbits, whose indicator and bifurcation curves are shown 
in Fig. 21 for the parameter space region of interest. For 
approximately a < 0.92, the curves coincide, resulting in 



0.918 0.92 
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FIG. 4: Indicator curves (icl and ic3) and bifurcation curves 
(bcl and bc3) for 7/8 periodic orbits. 



the standard scenario fSec. IIIlUl . The lower curve, cor- 
responding to bc3 and ic3, is the threshold for collision 
(i.e., reconnection) of hyperbolic orbits, while the upper 
curve, corresponding to bcl and icl, is for the simultane- 
ous collision (annihilation) of elliptic and non-symmetric 
hyperbolic orbits. For a > 0.92 the curves separate, and 
around a = 0.925 the bifurcation curve bc3, crosses the 
indicator curve icl. 



Recall that 6 = $ 



7/8, i I 



is the bifurcation curve for 



orbits on the symmetry line Si, and b = ^'7/3^^(0) is the 
indicator curve. In the following we drop the subscript 
7/8 and the dependence on a. Also, the symmetry of the 
SNM implies that <I>i = <I>2 and $4 = $3. 

The separation gives rise to the following two se- 
quences: 

For 0.92 < a < 0.925, the above curves are in the 
following order: < $3 < vj/^ < $1. 

1. Whereas for small perturbations {b < 'is), only 
the up and down (outer) 7/8-chains exist (wmin > 
7/8), ijJmin reaches 7/8 and inner orbits are born at 
6 = *3 (Fig.Et). Initially, this inner chain has the 
dipole topology, with elliptic orbits on S3 and S4, 
and hyperbolic orbits not on any symmetry line. 

(Fig.Eb). 

2. Somewhere in the range ^'3 < 6 < $3, a recon- 
nection between these off-symmetry line hyperbolic 
orbits with the outer hyperbolic orbits on S3 and 
S4 occurs (Fig. |S|:). At the point of reconnection, 
^max is 7/8 and continues to be so through step 
121 After this reconnection, meanders are born, the 
inner and outer chains display a nested topology, 
and Lu„,in < 7/8 (Fig.EJl). 

3. The inner elliptic orbits collide with the outer hy- 
perbolic orbits on the S3 and S4 symmetry line at 
6 = <I>3 (Fig.|S^), leaving only the off-symmetry line 
inner hyperbolic orbits and the outer elliptic orbits 
on the symmetry lines si and S2 (Fig. |Sf). 
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FIG. 5: 7/8-orbit reconnection at a = 0.923 (discussion in text). The upper plot shows a y vs. a; phase space plot {x £ 
[—0.25, 0.25], J/ G [—0.5,0]) and the lower one the corresponding winding number profiles along the two indicated symmetry 
lines {io along Si vs. y, lj e [7/8 - 0.0007,7/8 + 0.0005]). The plots correspond to the b-values: a) 0.22325262 (« ^3), b) 
0.22335, c) 0.22346060, d) 0.22355, e) 0.22364169 ($3), f) 0.22385, g) 0.22416343 (^i), h) 0.22465. 



4. At b ~ 'fi, the inner hyperbolie orbits eoUide 
(Fig. El;) and move onto the symmetry lines si and 
S2 (Fig.EJi). 

5. At 6 = <i>i, there is an hyperboHc-elliptic coUision 
on Sl and S2, after which 7/8-periodic orbits cease 
to exist and Wmax < 7/8. 



For a > 0.925, the order of the curves is changed to 
^3 < "S/i < $3 < $1. Most steps of the above sequence 
remain the same (with replacing $3 in the second 
step), except for the following: 

3. At 6 = ^I^i, the inner hyperbolic orbits collide and 
move onto the symmetry lines si and S2, while the 
hyperbolic elliptic orbits on S3 and S4 merely con- 
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tinue to approach each other. 

4. The mner eUiptic orbits coUide with the outer hy- 
perboUc orbits on S3 and S4 at 6 = $3 , leaving only 
the hyperbolic and elliptic orbits on si and S2- 

The difference between the two scenarios is the order in 
which the collision along S3 and S4 and the move of the 
hyperbolic points onto si and S2 occur. Both occur si- 
multaneously for the (a, 6)-value at which $3 and in- 
tersect. 

As in the standard scenario of odd-period orbits, me- 
anders appear in this case when the maximum of the 
winding number profile is rational. But now there are 
two meandering regions, and none of the meanders are 
are C/-invariant (Fig.^t, along S4). 

Note that the indicator curves, which are obtained us- 
ing the symmetries of the map, track collisions occurring 
along the ^/-invariant curve while the bifurcation curves 
track the collisions occurring at maxima in the winding 
number profile, which may or may not be shearless. 

This motivates an extension of the definition of bi- 
furcation curves. Recall that in the original definition 
the bifurcation curve determines the global boundary in 
parameter space between existence and non-existence of 
periodic orbits of the corresponding winding number. 
These curves correspond to maxima in the winding num- 
ber plot, and will be denoted by from now on. 

In the region of parameter space where four (or more) 
orbits of a particular winding number exist, a new kind 
of bifurcation curve can be defined that corresponds to 
parameter values at which the inner orbits are born, de- 
noted by In this case, the winding number profile 
has a local minimum (or local maximum for higher order 
curves). For even-period orbits, as expected and verified 
numerically, — . To estimate the region of pa- 
rameter space in which there are more than two chains 
of periodic orbits, we plot the points of separation of in- 
dicator and bifurcation curves for several different (odd 
and even) winding numbers. The results are shown in 
Fig. El Within numerical uncertainty, close to the 1/1- 
bifurcation curve the points seem to lie on a curve. Sim- 
ilar "curves" have been numerically observed for a few 
other winding numbers, e.g., Fig.[71 where we mark the 
points at which the indicator curves split from the bi- 
furcation curves (along si) for several inner orbits, in the 
region of parameter space close to the 1/3-orbit (left) and 
2/3-orbit (right) bifurcation curve (also shown). 

Remark 7. In Ref. Q Petrisor et al. conjectured that 
in generic nontwist maps the dipole formation scenario 
for even-period orbits does not occur. Instead, their nu- 
merical studies showed that one of the elliptic chains bi- 
furcates, creating and subsequently destroying a saddle- 
center pair, which has the effect of aligning the elliptic 
points of one chain with the hyperbolic points of the other 
one. Reconnection occurs then according to the standard 
scenario for odd-period reconnection. 



B. Non-standard reconnection for odd-period 
orbits 

Since there is no analog of the criterion Eq. (jSJ for 
the odd-period orbits, we cannot define indicator curves 
for these orbits. But we note that the indicator curves 
for even-period orbits are the bifurcation curves for the 
inner orbits, ^'i = Thus, we have implemented a 

version of the numerical method for finding bifurcation 
curves 0| for both the inner and outer orbits in the 
parameter space region where more than two chains exist. 
Fig. shows the bifurcation curves for the outer orbits, 
b = ^"s/g i(a)- and inner orbits, b = ^'^^g j^{a), of winding 
number 8/9. (Recall that in the odd-period standard 
scenario the bifurcation curves, are the same for all 

four symmetry lines.) For the remainder of this section, 
we will drop the dependence on a and the subscript 8/9 
for brevity. 

A magnification of (a,6)-spacc around the separation 
point (Fig. IHl reveals the birth of a new hyperbolic- 
elliptic pair of periodic orbits along each of the sym- 
metry lines at = <I>™ for a > 0.9295 (lowest curve in 
Fig. IHJ- This pair of orbits moves apart and eventually 
collides/annihilates with the outer, up and down, orbits 
on the symmetry lines. These collisions do not occur 
simultaneously in parameter space (in contrast to the 
even case), which explains the existence of two outer bi- 
furcation curves along each symmetry line, one for the 
outer up orbit and one for the outer down orbit, de- 
noted by (l)°"*'"P and ^o^'-'^o™^ respectively. Because 
of the symmetry of the SNM, <S>°''^;,^a^2 = ^s^andT" and 
and 2 ' = *3 and''4- To illustrate this non-standard 
reconnection-coUision sequence, we discuss in some de- 
tail the case of 8/9-periodic orbits along si at a = 0.94 
(Fig.EJ: 

1. Whereas for small perturbations {b < $'"), only the 
up and down (outer) orbits exist and LOmin > 8/9, 
inner chains of orbits are born (Fig.lHK) at 6 = 
and subsequently form a nested topology with me- 
andering orbits (Fig. ^p). As seen in the cor- 
responding winding number profile along si, the 
bottom of the valley exhibits a hill whose edges 
have the value Wmin = 8/9. The left edge, cor- 
responding to the new periodic orbit pair, broad- 
ens as its hyperbolic and elliptic orbits start to 
move apart (Fig. I10|l . whereas the right one, cor- 
responding to the heteroclinic connection between 
new pairs on other symmetry lines, continues to 
touch Wmin = 8/9 at only one point. The wind- 
ing number of the shearless curve, corresponding 
to the Zoca/ maximum (top of hill), and those of the 
meandering orbits are greater than 8/9. 

2. Somewhere in the range < 6 < (I>°"''"Pj the 
nested orbits reconnect (Fig. IHt) , resulting in four 
regular Poincare-Birkhoff chains (Fig. EJi), two on 
each side of the central shearless curve. 
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FIG. 6: Boundary of break-up diagram (dots), indicator curves (ic) and bifurcation curves (be) for two different winding 
numbers. Points marked (•) correspond to points where <^'^^ and $i" of close-by periodic orbits separate. 
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FIG. 7: Boundary of break-up diagram (dots) and bcl for 
1/3-orbits (top) and 2/3-orbits (bottom). Points marked (•) 
correspond to points where <I>5"' and $1" of close-by periodic 
orbits separate. 




FIG. 8: Magnification of split region in parameter space for 
8/9-bifurcation curves. The curves appear in the order < 

T,out,up ^ 3:out,down / i. . . j. j.\ 

3>i < "Pj^ . (discussion m text) 



3. At a higher 6-valuc in the same range, the two up 
chains and two down chains reconnect (Fig.|5^), re- 
sulting in two regions of nested orbits with mean- 
dering tori (Fig. 13^). Again, the winding numbers 
of the meanders is less than that of the reconnect- 
ing orbits (8/9). Each of these regions contains a 
shearless, but not ^-invariant, meander. 

4. At 6 = (I)°"*'"P^ the outer up orbits undergo a 
hyperbohc-elliptic collision (Fig. . 



5. At 6 = ([)°^*'down^ ^Y^^ outer down orbits undergo 
a hyperbolic-elliptic collision (Fig. Oi) , after which 
no further 8/9-periodic orbits exist. 
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FIG. 9: 8/9-orbit reconnection at a = 0.94 (discussion in text). The upper plot shows a, y vs. x phase space plot (x £ 
[—0.25, 0.25], y G [—0.5,0]) and the lower one the corresponding winding number profiles along the two indicated symmetry 
lines {u! along Si vs. y, lj € [8/9-0.0020, 8/9 + 0.0016]). The plots correspond to the 6-values: a) 0.2015781 (f=s $'"), b) 0.2018, 
c) 0.202155, d) 0.2036, e) 0.20515, f) 0.2059, g) 0.2067917 $°"''"P), h) 0.2081084 (f» $°"t.d°wn^j 



IV. DISCUSSION 

To understand the effect of the different reconnection- 
collision scenarios on the global transport properties of 
the map, we compare the computed bifurcation and indi- 
cator curves with the break-up diagram, obtained using 
the procedure of Sec. Ill El The results are presented in 



Fig.E] We see that these curves serve as a scaffoldiiig for 
the diagram, an observation also made in Refs. [ssUll ]. 
We will illustrate this interplay between the different re- 
connection scenarios and transition to global chaos by 
looking at the magnification of the break-up diagram 
around the indicator and bifurcation curves of the 5/6- 
orbits (Fig. 112). 
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FIG. 10: Magnification of the winding number profiles along 
si of Figs. |2t (top) andlHJj (bottom).. 



• For approximately a < 0.883, indicator and bifur- 
cation curves along the different symmetry lines co- 
incide, and lie within the non-chaotic region of the 
break-up diagram. Fixing a and increasing the per- 
turbation b through the rcconncction process does 
not lead to global chaos. This is because of the 
presence of invariant tori "above" and "below" the 
reconnecting chains of orbits. (For the remaining 
discussion, we will refer to these as "twist" tori.) 

• For approximately 0.883 < a < 0.893, the lower 
pair of curves follows the boundary of the break-up 
diagram. Increasing the perturbation through the 
boundary corresponds to a rcconncction that re- 
sults in global chaos. The reason for this can be un- 
derstood as follows: For parameter values slightly 
below the lower pair of curves, no twist invariant 
tori exist above or below the shcarless region. Only 
the tori in the region between the two chains of 5/6- 
orbits inhibit transport. When the perturbation b 
reaches the lower pair of curves, the hyperbolic 5/6- 
orbits collide and the tori between the two chains 
are destroyed. Since the colliding orbits are hyper- 
bolic, their separatrices form a heteroclinic tangle 
that connects the up and down chaotic regions re- 
sulting in global chaos, consistent with the observa- 
tion that the lower pair of curves form a boundary 
of the break-up diagram. 

• For approximately 0.893 < a < 0.91, we see the 
first of the two non-standard scenarios (Fig.lSJ dis- 



cussed in Sec. IIII Al Here again, no twist tori exist 
for parameter values close to the lowest curve 4*3. 
The inner orbits, which are born at b — '^3, move 
apart and reconnect with the outer hyperbolic or- 
bits for some ^3 < 6 < $3. At this point, global 
chaos ensues because the up and down chaotic re- 
gions are joined through the reconnecting invariant 
manifolds. Thus the lower boundary of the break- 
up diagram lies between \E'3 and $3 . 
Also, for 5 > Vfi, the off-symmetry line hyper- 
bolic orbits move onto the si symmetry line and 
their invariant manifolds are no longer connected 
("de-connection"). Meandering orbits are created 
as these hyperbolic orbits move apart along si. 
These meandering orbits (Fig. inhibit global 
chaos. Thus, 6 = (seen as the threshold for 
"de-connection" ) forms the upper boundary of the 
break-up diagram. 

• It was noted in Ref. [s^ that if the up and down 
twist regions of the map exhibit chaos when the 
odd-period orbits reconnect, the rcconncction leads 
to global chaos for the above-mentioned reason 
(i.e., because of the heteroclinic tangle of the in- 
variant manifolds). 

Thus we see that, in general, reconnection (and the 
break-up of twist tori) has a significant effect on the tran- 
sition to global chaos. Many, if not all, of the "smooth" 
boundaries of the break-up diagram can be conjectured 
to be the result of reconnections of invariant manifolds. 
Indeed, the reconnections of orbits with higher periods 
can give rise to finer structures in the break-up diagram. 
Thus, the "fractal" nature of this diagram might be re- 
lated to the reconnections in phase space happening at 
smaller and smaller scales. A reliable criterion for de- 
termining the reconnection threshold and its numerical 
implementation will be extremely useful for testing these 
ideas. 



V. CONCLUSION 

In this paper we presented a unified view of various 
bifurcations and reconnections that occur in the stan- 
dard nontwist map, but are not observed in twist maps. 
The bifurcations are locally of three types: collision and 
annihilation of hyperbolic-elliptic orbit pairs, collision of 
two symmetric hyperbolic orbits resulting in two non- 
symmetric hyperbolic orbits, and the simultaneous (in 
parameter and phase space) collision and annihilation of 
two hyperbolic and two elliptic orbits. The latter is seen 
only in the SNM (and closely related maps) because of 
its high degree of symmetry. The two former bifurcations 
are observed in all nontwist maps. 

We have demonstrated that the non-generic standard 
nontwist map has more types of periodic orbit recon- 
nection scenarios than previously known, resulting from 
the presence of four (or more) chains of periodic orbits 
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FIG. 11: Parameter space showing the points for which shearless invariant tori exist. Also shown are some bifurcation and 
indicator curves and the critical points (marked by o) found using Greene's residue criterion. 
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FIG. 12: Section of parameter space around the 5/6-orbit indicator and bifurcation curves. 



of the same winding number. Meandering tori, associ- 
ated with the odd-period reconnection in the SNM, have 
been shown to exist for certain regions of parameter space 
for non-standard even-period reconnection. Petrisor and 
co-workers have started to compile a list of possible 



reconnection scenarios for nontwist maps. Our investi- 
gation contributes new global scenarios for the case of 
nontwist maps with a reversing symmetry group, includ- 
ing a spatial symmetry. 

We also conjecture (and present heuristic observations) 
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that the reconnection of hyperbolic separatrices leads to 
global chaos if no "twist" tori exist in the region "out- 
side" the reconnecting chains of orbits. Determination of 
the break-up threshold of the last twist tori and a pre- 
cise criterion for calculating reconnection thresholds will 
shed more light on the accuracy and usefulness of these 
observations. 
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APPENDIX A: BASIC DEFINITIONS 

For reference, we list a few basic definitions used 
throughout the main text. An orbit of an area-preserving 
map M is a sequence of points {{xi, ?/i)}j-_oo such that 
M {xi,yi) = (xi+i, j/i+i). The winding number lo of an 
orbit is defined as the limit lo = Imii—yaoixi/i)^ when it 
exists. Here the x-coordinate is "lifted" from T to M. 
A periodic orbit of period n is an orbit M" (xi,yi) = 
[xi +m,yi), Vi, where m is an integer. Periodic orbits 
have rational winding numbers lo = m/n. An invariant 
torus is a one-dimensional set C that is invariant under 
the map C = M{C). Of particular importance arc the 
invariant tori that are homcomorphic to a circle and wind 
around the .T-domain because, in two-dimensional maps, 
they act as transport barriers. Orbits belonging to such 
a torus generically have irrational winding number. 

APPENDIX B: SYMMETRY PROPERTIES OF 
THE SNM 

In this appendix, we review the symmetry properties 
of the SNM. A detailed discussion of these types of sym- 



metries in general can be found in Ref. [43 and, in the 
context of nontwist systems, in Ref. 0. 

The standard nontwist map M is time-reversal sym- 
metric with respect to two involutions Iq and /i, i.e., 
Af-i = /ri M If ^ Id, for i = 0,1. It then fol- 
lows that the SNM can be decomposed as M ~ Ii o Iq. 
In addition, Al is also symmetric under an involution 
S which commutes with both li, i.e., M = S^^ AI S, 
S'^ = Id, and S li = li S. The involutions have the fol- 
lowing form: /o(a;,j/) — {~x,y — &sin(27ra;)), Ii{x,y) ~ 
{~x + a(l - y^),y), and S{x, y) ^ {x + 1/2, -y). 

Since the li arc orientation reversing, their fixed point 
sets, Ti = {{x,y)\Ii{x,y) = (a;,y)}, arc one-dimensional 
sets, called the symmetry lines of the map. For the SNM, 
To consists of si = {{x,y)\x = 0}, S2 = {{x,y)\x = 1/2}, 
while Fi consists of S3 = \^{x,y)\x = a (l — y"^) /2}, 54 = 
{{x,y)\x = a{l-y^)/2+l/2}. 

It is easy to check that S does not have any fixed 
points. Symmetry lines are useful because the numeri- 
cal search for symmetric periodic orbits, i.e., orbits with 
a point belonging to one of the symmetry lines, is a one- 
dimensional root finding problem and hence consider ably 
easier than the search for non-symmetric orbits. [i^ l5Cl|| 
For the SNM, there are generally two orbits along any 
symmetry line of any winding number, called the up and 
down orbits. 

The SNM is time-reversal symmetric with respect to 
maps Sli as well. The fixed point sets of these orientation 
preserving maps consist of points in phase space, called 
indicator points. For the SNM, fixed points of 5'/o and 
SIi are, respectively, z'^'^^ — (Tj^T^), and z^^'^^ = 
(fTi,0). 

The involutions /q, /i, S, and their products form a 
group Q called the reversing symmetry group. There is at 
most one homotopically nontrivial invariant torus that is 
invariant under the action of This torus is called the 
central or Q -invariant shearless torus 7^. The indicator 
points belong to 75 when it exists. 
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